Signature of the interaction between dark energy and dark matter in 

observations 

Elcio AbdallaQ L. Raul Abramo[|] and Jose C. C. de Souz^^l 

Instituto de Fisica, Universidade de Sao Paulo, 
CP. 66318, 05315-970, Sdo Paulo, SP, Brazil 

Abstract 

We investigate the effect of an interaction between dark energy and dark matter upon the dynamics 
of galaxy clusters. This effect is computed through the Layser-Irvine equation, which describes how an 
astrophysical system reaches virial equilibrium and was modified to include the dark interactions. Using 
observational data from almost 100 purportedly relaxed galaxy clusters we put constraints on the strength 
of the couplings in the dark sector. We compare our results with those from other observations and find 
that a positive (in the sense of energy flow from dark energy to dark matter) non vanishing interaction is 
O , consistent with the data within several standard deviations. 

cr 

tlj PACS numbers: 98.80. C9; 98.80.-k 

OJQ, 

m 

> 

CO 
CN 

in 
d 

On 
O 

> 

X 



'Electronic address: eabdalla@fma.if.usp.br 
' Electronic address: abramo@fm a.if. usp.br | 
^Electronic address: jccampos@fma.if.usp.br 



I. INTRODUCTION 



One of the most surprising results of physics and cosmology in the last ten years is the cos- 



mological accelerated expansion, which has been proved beyond reasonable doubt by observations 
[l-^]. The simplest explanation of such observations is a cosmological constant, which is, however, 
off theoretical computations by 120 orders of magnitude, a result that calls for urgent explanations. 
Another possibility is to argue that the Universe contains a strange dynamical component with 
negative pressure, dark energy, which is responsible for more than three quarters of its energy 
content 4] . 

Interactions of dark energy or dark matter with baryonic matter and radiation must be either 
inexistent or negligible in order to comply with stringent observations of the visible sector. However, 
in a field-theoretical approach, dark matter is some particle of a unification scheme, and dark 
energy should also be described by a field, which means that some level of interaction between 
these different sectors is basically mandatory. In the literature there is a large body of work 
dealing with such a possibility, e.g., It is interesting to notice that a coupling between dark 



energy and dark matter can also serve to alleviate the coincidence problem 



111 ]. The value of 



the coupling and holographic arguments also allow for the crossing of the phantom barrier which 



separates models with equations of state w = pip > —1 from models with w < —1 [131 _ see 

nn 

also pJ,Ha|. Further consequences of the interaction have been studied, as e.g. its effect on the 



lowest multipoles of the cosmic microwave background (CMB) spectrum 



13, 
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161 ]. The strength 



. Comparison with 



of the coupling is presumably as small as the fine structure constant 
supernova data and CMB and large-scale structure have also been analysed [18(. Nevertheless, the 
observational limits on the strength of such an interaction remain weak [lj]] . 

In addition to these constraints, dynamical dark energy (i.e., a time- and space-dependent field) 
has an impact on the large-scale structure due to its fluctuations. In that case, dark energy affects 
the pro cess of structure formation by means of its density fluctuations, both in the linear 
231 and the non-linear 
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affected 



3, 
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26] 
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11}, 



251 ] regimes, as the growth of dark matter perturbations can be 



However, most constraints on dark interactions concern the asymptotic behaviour of the dark 
sector, both in time and in space. Local or nearby checks are still weak. This started to change 
with the detailed analysis of galaxy clusters and their internal structure. It was suggested that 
the dynamical equilibrium of collapsed structures could be affected by the coupling of dark energy 
to dark matter, which as a result could affect the averaged energy distributions predicted by the 
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virial theorem - in fact, by its relativistic counterpart, the Layser- Irvine equation [27(. This was 



281] . The basic idea is that the virial 



first proposed, and analyzed for the Abell cluster A586, by 
theorem is distorted by the mass non-conservation generated by the coupling of dark matter with 
dark energy. 

□ 

In a previous paper 29] we showed how the Layser-Irvine equation, describing the evolution 
towards virialization, is changed by the presence of the coupling. We showed that this violation 
leads to a systematic bias in the estimation of virial masses of clusters when the usual virial 
conditions are employed. By comparing weak lensing and X-ray mass-observables on the one hand, 
and virial masses on the other hand, we were able to cross-check not only the strength of the 
coupling, but also whether it was indeed a feature of the virial mass estimate, and not a more 
complicated set of biases between all these mass-observable relations. Our main result was that 
a single parameter could explain all the bias between virial mass and the other estimates (weak 
lensing and X-ray) to 95% confidence level (C.L.), or, equivalently, two standard deviations. 

However, the amount of data (about thirty independent useful observations) did not allow a 
detailed numerical analysis, resulting in relatively weak constraints on the coupling. In this paper 
we will not only reanalyze that data, but will also include a new set of data for almost a hundred 
X-ray observations of galaxy clusters, and compare the presumed mass from those observations 
with the virial masses for those clusters. 

Notice that even though the uncertainties associated with any individual galaxy cluster are 
very large, by comparing the naive virial masses of a large sample of clusters with their masses 
estimated by X-ray and by weak lensing data, we are able to constrain the bias between them and 
to impose much tighter limits on the strength of the coupling than has been achieved before. 



II. PHENOMENOLOGY OF COUPLED DARK ENERGY AND DARK MATTER MOD- 
ELS 

We start with a simplified two fluid model interaction: 

p dm + 3Hp dm = ip (1) 
p de + 3Hp de (l + w de ) = -ip , 

where a dot denotes time derivative, H is the expansion rate, p dm and p de are respectively the 
energy densities of dark matter and dark energy, and w dm and w de are their equations of state. 
Notice that the continuity equation still holds for the total energy density pTot = Pdm + Pde- 
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Phenomenologically, one can describe the interaction between the two fluids as an exchange of 
energy at a rate proportional to a combination of the energy densities 

<A = CHpTot ■ (2) 
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311 ] . but that does not concern 



In fact, the term proportional to p dm could lead to instabilities 
us here in view of the local character of the computation. We are interested in collapsed structures 
where the local, inhomogeneous density a is far from the average, homogeneous density p. In that 
case the continuity equation for dark matter reads 251 ]: 



<?dm + 3Ha dm + V {a dm v dm ) = (H (a dm + a de ) , (3) 

where v dm is the peculiar velocity of dark matter particles. We have considered the local density 
of dark energy to be proportional to the local density of dark matter, a de = b em a dm . If for a given 
model the dark energy component is very homogeneous, then b em ~ 0. We do not consider the case 
where b em depends on the size and mass of the collapsed structure - although this should probably 
happen in realistic models of structure formation with dark energy [25J. Hence, the continuity 
equation with dark matter coupled to dark energy reads: 

& dm + 3Ha dm + V (a dm v dm ) = (Ha dm , (4) 

where C = C(l + b em )- 

Considering that the kinetic (K dm ) and potential (U dm ) energies of a set of particles interact 
via gravity and the coupling ([2]), we found that the corrected Layser- Irvine equation leads to the 
condition: 

(2-()K dm + (l + b em )(l-2()U dm = . (5) 
Taking £ = b em = we recover the usual virial condition. 

III. MASS ESTIMATION AND LIMITS ON THE COUPLING 



In order to find deviations from the virial t 
largest virialized structures in the Universe 32J, 



reorem we analyse galaxy clusters, which are the 



331 ] . Studies of cluster in the optical, in X-ray and 



through weak lensing are available in the literature, which can then be used to estimate the mass 
of the cluster, under some hypotheses. 

Clusters can be observed directly in the optical wavelengths, where the radial velocities are 
estimated through their projected (line-of-sight) velocities. The velocity field, together with the 



projected spatial distribution of galaxies in the cluster, allow for an estimate of the relative shares 
of the potential and kinetic energy of its constituents. Assuming that the clusters are virialized, 
their masses can then be computed. But if there is a coupling of dark energy to dark matter, we 
find that [29| 



(l + b em )^ = -2 1 -4B. (6) 
K dm 1 - ZQ 

Since the potential energy is proportional to the square of the masses, and the kinetic energy is 
only linearly proportional to the mass, there is a shift in the usual virial mass estimation, which is 
entirely due to the coupling, by a factor of (1 — 2£)/(l — C/2)- 

One can measure the mass with other observations by making very simple assumptions which 
have nothing to do with the precise nature of the equilibrium of the system. In particular, both 
the weak lensing and X-ray methods provide physical observables which can be used to estimate 
the total mass of a cluster in ways which are totally independent of the virial method. 

In what follows we will assume that the variations in the mass estimators can be explained 
by a single parameter - the coupling between dark energy and dark matter. Of course, all these 
mass estimates are rife with systematic errors of many types, and intrinsic uncertainties in these 
estimates should be expected. For instance, if the gas in a cluster is not homogeneously distributed, 
the clumps of gas will enhance the X-ray luminosity of the cluster relative to a homogeneous 
distribution (as the X-ray luminosity is proportional to the square of the number density of free 
electrons), leading to a possible overestimation of the X-ray mass. Clearly, our method cannot make 
a distinction between the constraint on the coupling of the dark sector and some other source of 
systematic error, such as the dumpiness of gas in clusters. 

We parametrize this single-parameter deviation in mass estimations, both with weak lensing or 
X-ray data, as: 

M mr = ^S^M X = ^B MwL (7) 



Previously 29J, we compared and cross-checked the three ratios that can be obtained from each 
pair of masses. We obtained that £ = 0.04 ± 0.02, and checked that there was no bias between 
X-ray and weak lensing estimates. 



The data in 34] allows to more than triple our dataset, compared with what we had before, 
leading to the possibility of a more thorough numerical analysis. First, we can achieve a better 
accuracy, showing that the coupling is non zero by several standard deviations. And second, we 
also analyzed a possible dependence of the coupling on the size of the cluster as well as on its 



redshift - checking for possible evolution effects. Hence, since many more X-ray data for clusters 
are available than for weak lensing, we can enhance our test of the dark interaction by including 
X-ray clusters without the weak lensing counterpart. 

We can collect all these datasets by taking each mass estimate, together with its uncertainties, 

iinction. Our data comprises a set of 34 weak lensing, X-ray and 



and constructing a likelihood 
virial mass estimates from as well as 61 X-ray and virial (optical) mass estimates from 

34|. The weak lensing and X-ray masses were compared with the optical data, and the results are 
represented in figure ([T]). 

One of the difficulties, already encountered in our previous paper, is how to deal with asymmetric 
errors. At that time we just worked with the average error, discarding corrections due to such an 
asymmetry. Taking this asymmetry into account, e.g. by making use of the techniques presented 
in Ref. jfj], may lead to more accurate constraints, however, there are some difficulties associated 
with the asymptotic behaviour of the likelihood function, and the final results are not significantly 
different from the ones obtained with a more simple analysis. Hence, we employ the simpler, 
symmetric treatment of the errors. 

With the 61 X-ray masses mentioned above (including, of course, the optical data), as well as 
the earlier 34 X-ray and weak lensing masses, we can construct the likelihood function C for Q with 
the full dataset: 



N 

C oc exp 



1 



2 of 



1-2C , 

Ji 



l-C/2 



(8) 



i=l 

where / represents the mass ratios Mx/M v i r and Mwt/M v i r , and the errors a are approximated 
by a simple geometric mean of the asymmetric uncertainties for each of the data, ^/cr + (T_. Setting a 
conservative prior for the values of £ to be flat, and between -0.3 and 0.3, we obtain the probability 
density function displayed in figure ((TJ). We learn that X-ray data are more robust. The weak 



lensing data are more spread, though consistent with X-ray measurements [291 ] . 
The joint analysis point to: 

C = 0.14 ±0.01 . (9) 

Notice that the likelihood function assumes implicitly that all the variation in the mass esti- 
mations can be attributed to a single parameter, and the peak and width of that likelihood give 
us the best-fit and the uncertainty for that parameter. We could also have made no assumption 
whatsoever about the nature of these estimates, and compute the sample variance of all mass esti- 
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Figure 1: Probability density function (p.d.f.) for ( obtained with all available cluster masses. The result 
points to ( = 0.14 ± 0.01, that is, 14 standard deviations from the null result. 



mates. The sample variance of the ratio between mass estimates (which is, again, degenerate with 
our coupling parameter) results, of course, in a much higher variance, namely a 1 = 0.3. 

The main result obtained here is that, assuming that a single parameter can explain the dif- 
ference between the mass estimations, there is a positive coupling which is nonzero within several 
standard deviations. Now, in order to investigate whether different masses or evolution effects, for 
example, could be playing a role in our analysis, we can divide our data into subsets. 

First, we can redo the numerical analysis after dividing the data in terms of the velocity disper- 
sion - which is in fact a proxy for mass. This analysis leads to the three diagrams shown in figure 




Figure 2: P.d.f s concerning the cluster data divided in three sets, in increasing order of the dispersion 
velocity (up to 600 km/s, from 600 to 1000 km/s, and above 1000 km/s.) The data are compatible with one 
another and the results point to essentially the same values of all at least three standard deviations from 
the null result. 



7 



The results for £ are, respectively, for smaller to larger values of the velocity dispersion: 

( Vl = 0.14 ±0.04 

( V2 = 0.17 ±0.01 (10) 
C„ 3 = 0.13 ±0.01 . 

We can also check how robust are these results concerning the redshift - that is, we could ask 
whether evolution effects are important. Hence, we divide the data in subsets in order of increasing 
redshift - and we do this only for the X-ray observations. The results are displayed in ©. 





Figure 3: P.d.f's concerning the cluster data divided in two sets, in increasing order of the redshift (below 
and above z = 0.05). The data are compatible with one another and the results point to compatible values 
of £, at several standard deviations from a vanishing result. However, in the absence of weak lensing results 
the results point to a considerably higher value of the coupling compared with the previous ones. 



The results for ( are, respectively, for smaller to larger values of the redshift: 

( Zl = 0.24 ±0.01 

Cz 2 = 0.34 ±0.01 . (11) 

Finally, some words are in order about choices made for the available X-ray data for cluster. Not 
all clusters displayed in 



34| are relaxed. Moreover, some data have different results for different 
kinds of observations. In order to deal with these differences we used two procedures. In one, we 
always considered the data with largest radii. Since it is just a matter of random choice, the actual 
result, if robust, should not depend on it. In the second criterium, we considered the smallest radii. 
The results obtained using smaller and larger values of the radius are, respectively: 

C n = 0.32 ±0.01 

Cr 2 = 0.30 ±0.01 . (12) 

Hence, we conclude that for both criteria the results of our analysis are very similar, which means 
that our constraints are sufficiently robust with respect to these choices. 
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IV. COMPARISON WITH PREVIOUS RESULTS 



In addition to the constraints obtained above, there are several other results concerning the 
interaction of dark energy and dark matter. A possible transition of the Dark Energy equation of 
state has been observed and an analysis of an explanation of such a transition in terms of Dark 
Energy and Dark Matter interaction in the terms proposed here has been performed in Such 
an analysis leads to results which can be translated to constraints in £ as 

C = 0.18 ±0.18 . (13) 

That study has been made using a holographic model for Dark Energy, but essential was the 
interaction between the two dark sectors. 

Furthermore, a study of the age of the old quasar APM 0879 + 5255 also constrains the inter- 
action Indeed, such a quasar is too old to be accomodated by the standard cosmology, and 



we have shown in 



131 ] that an interaction leaking energy from the dark energy into dark matter 



sector such that the coupling constant is given by 

C = 0.36 ±0.18 (14) 

naturally accomodates the age of the quasar being also compatible with the age of the universe as 
given by CMB observations. 



381 ] deviation from the 



On the other hand, the small £ behavior of the CMB angular spectrum 
standard model leads to further constrains in the interaction as given by [lj] 

C = 0.45 ±0.15 . (15) 

Furthermore, the interaction can also be modeled by an alternative interacting field theory 
lagrangian, where Dark Energy is described by a tachyonic lagrangian while the role of dark matter 
is played by a fermionic field. A minimal interaction of the Yukawa type completes the picture. 
Such a model is constrained by baryon acoustic oscillations (BAO), lookback time, supernovae and 



CMB shift parameter 



. The effective coupling of both sectors, namely fermionic (dark matter) 



and bosonic (dark energy) is constrained to fulfill 

C = 0.17 ±0.09 . (16) 

Further observational data (once again taking into account supernovae, shift parameter, BAO 
and galaxies ages estimates data) together with a phenomenological model based on a two fluid 
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interaction also points to some (feeble) constraints 40J] 



C = 0.01 ±0.04 , (17) 

where we have considered only the most conservative result of that analysis (in fact, the results 
present a weak dependence on the type of interaction; we have taken the average, most conservative 
result). 

Entirely theoretical arguments derived from a thermodynamical analysis allows an estimation 
of the interaction parameter for the two fluid model as 41], 



C = 0.15 ±0.15 . (18) 

We also recall here a previous result concerning a smaller set of clusters (contained in the total 
data set used in the present work): {29]. 

C = 0.04 ±0.02 . (19) 

Our joint set of cluster observations led us to the constraint of Eq. Q. This is one way of viewing 
our result. However, considering the subsequent analysis (in particular the sample variance), it 
may be fair to see the errors as much larger, possibly of the order of magnitude of the highest 
difference in sub-averages. Considering the highest spread to be 5( = 0.5 x (0.34 — 0.14) = 0.10, 
we take our results here as: 

C = 0.14 ±0.10 . (20) 

Even in this pessimistic case, we can infer that our result supports the previous ones, i.e., that 
the interaction coupling parameter is a small but positive and of the order of 0.1, and almost three 
standard deviations from vanishing. 



V. CONCLUSIONS 



We are fully aware that our lack of precise knowledge about the errors arising from systematics 
could weaken our conclusions. We recall here that our main work hypothesis is that the differences 
between the various kinds of masses estimates come solely from the fact that they do or do not 
take into account the interaction parameter. Obviously, this assumption is not absolutely true, 
but represents an approximation we use to replace the need to analyze any particular systematic 
effects. Also, we have used in this work the best available data sets well-suited for our purposes. 
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Nevertheless, we still believe that our analysis, combined with previous results, signals in a coherent 
manner towards the possibility of the dark energy-dark matter interaction. 

The scenario in which dark energy and dark matter are in interaction seems to be a strong 
physical possibility. The interaction intensity is a small constant, possibly larger than the fine 
structure constant, in agreement with [171 ] . The results are all consistent among themselves, and we 
obtain the interesting conclusion that the interaction constant is far from zero by several (possibly 
three) standard deviations, which is a rather intriguing conclusion. The positive value is also in 



421 ]. strengthening our conclusions. A few 



agreement with the thermodynamical arguments of 
words should be said about the possible influence of the interaction parameter on the cosmic 
history (in particular, in the structure formation process). The behavior of dark matter density as 
a function of time can be rather different from the non interacting case depending on the kind of 
interaction considered. We do not think that the proposed interaction should be the same along the 
whole history since decoupling, when a more robust suggestion has to be made. Since usually the 
matter-radiation equality is supposed to occur, in the standard model (e.g. in the non dark sector 
case) as well, the order of magnitude of the cosmological time parameter should not get strongly 
modified. (Notice that in the past, with interaction, dark matter density was smaller compared to 
the noninteracting case. Thus, the matter density is bounded from below by the baryonic matter 
density.) 
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